Draft version February 28, 2013 

Preprint typeset using I^TgX style emulateapj v. 5/2/11 



DANCING IN THE DARK: NEW BROWN DWARF BINARIES FROM KERNEL PHASE INTERFEROMETRYi 

Benjamin Pope 

Sydney Institute for Astronomy, School of Physics, University of Sydney, NSW 2226, Australia 



o 

tin 



6 

GO 



Frantz Martinache 
National Astronomical Observatory of Japan, Subaru Telescope, Hilo, HI 96720, USA 

Peter Tuthill 

Sydney Institute for Astronomy, School of Physics, University of Sydney, NSW 2226, Australia 

Draft version February 28, 2013 

ABSTRACT 

This paper revisits a sample of ultracool dwarfs in the Solar neighborhood previously observed with 
the Hubble Space Telescope's NICMOS NICl instrument. We have applied a novel high angular 
resolution data analysis technique based on the extraction and fitting of kernel phases to archival 
data. This was found to deliver a dramatic improvement over earlier analysis methods, permitting a 
search for companions down to projected separations of ~1 AU on NICl snapshot images. We reveal 
five new close binary candidates and present revised astrometry on previously-known binaries, all of 
which were recovered with the technique. The new candidate binaries have sufficiently close separation 
to determine dynamical masses in a short-term observing campaign. We also present four marginal 
detections of objects which may be very close binaries or high contrast companions. Including only 
confident detections within 19 parsecs, we report a binary fraction of at least = 17.2^^'^%. The 
results reported here provide new insights into the population of nearby ultracool binaries, while also 
offering an incisive case study of the beneffts conferred by the kernel phase approach in the recovery 
of companions within a few resolution elements of the PSF core. 

Subject headings: techniques: interferometric — techniques: image processing — techniques: high 
angular resolution — stars: low-mass — stars: formation — (stars:) brown dwarfs 
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1. INTRODUCTION 

A detailed picture of multiplicity is an essential ele- 
ment to understanding low-mass stars and brown dwarfs. 
Binary systems present an opportunity to determine 
model-independent dynamical masses when both astrom- 
etry and radial velocity data are available. Systems so 
characterized may then become part of the foundations 
for the construction of an observationally constrained 
mass-luminosity-age sequence applicable across the en- 
tire class. 

Furthermore, the statistical properties of populations 
of low mass binaries have profound implications on the 
basic physics of star formation and solar system assem- 
bly. Multiplicity rates are a key discriminant between hy- 
potheses about the formation and evolution of lo w mass 
systems, as discussed in Burgasse^^^2 (!2QQ7l ). Two 
main mechanisms have been proposed for the formation 
of brown dwarfs in t he field: ern bryo ejection, and gravo- 
turbulent collapse (jBasul l2Q12f ). Specifically, the em- 
bryo ejection hypothesis predicts a low binarity incidence 
(~ 8%) (Bate 2012) . which confiicts with the observed bi- 
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^ Based on observations performed with the NASA/ESA Hub- 
ble Space Telescope. The Hubble observations are associated 
with proposal ID 10143 and 10879 and were obtained at the 
Space Telescope Science Institute, which is operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under 
NASA contract NAS 5-26555. 



narity rate (- 15%) (iReid fc Hawlevll2QQ5[ ). Mapping the 
incidence of binarity, and in particular extending com- 
pleteness to smaller orbital separation is therefore of in- 
terest in establishing the primary formation mechanism 
of field brown dwarfs. 

Snapshot imaging is a straightforward way to discover 
new multiple systems. Intrinsically faint and red, L- 
dwarfs present challenging targets for ground based ob- 
servations, typically requiring laser guide star (LGS) 
adaptive optics (AO). Space telescopes naturally offer 
high Strehl ratio imaging at diffraction-limited resolu- 
tion, with the major departures from ideal performance 
arising from field-dependent PSF changes, spacecraft jit- 
ter, and slow optical drift from thermally induced breath- 
ing modes of the mechanical structures. 

Imaging campaigns with the HST have demonstrated 
notable success in prospecting for companions to cool 
objects, providing high quality diffract ion- limited images 
of a large number of targets (Reid et all l2QQ6al l2QQ8f ). 
These campaigns have shed light on the population of 
cool dwarfs in the Solar neighborhood. 

The simplest and most widely used method for detec- 
tion of companions in snapshot imaging essentially relies 
on direct visual examination of images. Obvious compan- 
ions are quickly identified, and traditional astronomical 
image analysis tools, namely aperture photometry and 
centroiding, provide the important astrometric and pho- 
tometric characteristics of the target. Faint or close-in 
companions are, however, easily missed in a visual search 
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and identifying such objects requires more sophisticated 
computational techniques. 

For example, some stell ar images exhibit an elongation 
along one axis as noted byl Reid et alJ (j2QQ6af ) which may 
be suggestive of the presence of a barely-resolved com- 
panion. Subtraction of a model PSF h as been exploited 
to infer the presen c e of a companion (|Krist et al.l 119981 : 
iPravdo et all I2QQ4I : iDieterich etld] l2Q12f ). although the 
performance of this approach is arguably poor, and fur- 
thermore it weakly constrains the relative photometry 
and astrometry. We propose to look at the same im- 
ages from an interferometric standpoint, leveraging the 
exquisite level of calibration this technique offers. 

For the detection and characterization of com- 
panions at small angular separations, non-redundant 
masking (NRM) interferometry used in conjunction 
with AO h as dem onstrated outstand ing perfo rmance, 
e.g. in iTuthih et a l. (2006); Lloyd et al.l (12006); 
iKraus & Ireland! (|2012f ) . The key underpinning such suc- 
cesses has been the robust, self-calibrating nature of the 
observables recovered from NRM interferometry, and in 
particular the closure phase^ first suggested for the ra- 
dio (Jennison 1958) and later exploited in the optical 
([Bald win et aJJ [19861). Imaging systems where the phase 
on any given baseline in the pupil is disturbed by ran- 
dom phase errors from atmospheric or instrumental aber- 
rations suffer from degraded performance. However by 
summing phases around closed loops of non-redundant 
baselines, these random phasors cancel out and the re- 
sulting closure phases are extremely robust to wave- 
front aberrations. N RM interferometry from the ground 
(jTuthill et al.l I2Q00I ) relies heavily on closure phase for 
high-contrast detection, and the re are plans to extend the 
technique to space platforms (Sivar amakrishnan et al.l 
l2009f ). Recent o bservations achiev e d with this tech- 
nique re ported bv iLlovd et al.l (120061): Martinache et al.l 
(I2007D: llreiand et all (120081 ): iKrauset al.. (2008,) and 
iMartinache et all ()2009f ) demonstrate that the level of 
calibration achieved with interferometric measurements 
permits the detection of companions at scales at or even 
somewhat beyond the diffraction limit of the imaging sys- 
tem. Recently, NRM interferometry succeeded in provid- 
ing evidence for a low- luminosity c ompanion in the tran- 
sitional disk h ost systems T Cha (jHuelamo et al.ll201ll ) 
and LkCa 15 (iKraus fc Ireland l2012[ r 

It has recently been demonstrated that if a conven- 
tional (full- aperture) PSF is of sufficient quality (wave- 
front residual errors typically < A/4), an analogous set of 
high-quality inter ferometric observa bles can be extracted 
from the images (jMartinachell2010f ). These new quanti- 
ties are the kernel phases^ and represent a generalization 
of the idea of closure phase to a redundant pupil config- 
uration. The major advance offered by kernel phase in- 
terferometry is that it is not restricted to non-redundant 
pupils. In brief, for small wavefront errors (i.e. high 
Strehl ratio), the phase errors in the pupil plane can be 
related to those in the Fourier plane by a linear opera- 
tor. The kernel or null-space of this operator therefore 
singles out a subspace of baseline phases which are not 
affected by this error, which can then play the same role 
as closure phases in providing a robust set of observ- 
ables to constrain image structure. Kernel phases were 
first successfully extracted from HST/NICMOS data on 
a single target by .Martinache (2010.) , demonstrating sig- 



nific ant improvement oy er more traditional data anal- 
ysis (jPravdo et al.ll2QQ3 ). The technique has now also 
been successfully applied to ground based AO observa- 
tions (Martinache 2011). 

This paper revisits a sample of nearby ultracool dwarfs 
observed by the Hubble Space Telescope NICMOS NICl 
camera and first presented in iReid et al.l ()2006a[ ) and 
iReid et aTl ()2008f ). Our analysis allows dramatic exten- 
sions to the discovery space for putative companions, and 
in particular explores separation ranges down to 1 AU on 
targets located within 20parsec. Section [2] provides an 
overview of the dataset and introduces the methods used 
for our new analysis. Section [3] discusses the results of 
the kernel phase analysis for the entire sample and im- 
plications for the astrophysical interpretation of brown 
dwarf formation. 

2. SAMPLE AND METHODS 
2.1. Sample of Ultracool Dwarfs 

This study focuses on two samples of ultracool dwarfs, 
observed with the HST/NICMOS NICl came r a, and 
who se properti e s wer e reported by iReid et aTl ()2006af ) 
and lReid et al.l (2008). Each target was observed in two 
filters: FllOW and F170M, which correspond loosely to 
the astronomical J and H bands. These differ in that 
the J and H bands sample atmospheric transmission 
windows, which do not constrain space-based observa- 
tions. We will use J and H as shorthands for FllOW 
and F170M respectively, but the difference should be 
noted. Table [T] summarizes the observational properties 
of the combined sample as stated in Tab le 1 of Reid et al.l 
([2006a ) and Table 1 of lR.eid et al.l (f2008>) . 

In addition to detecting several binaries by traditional 
data analysis methods, these authors also provide a list 
of 43 and 26 apparently unresolved objects in the 2006 
and 2008 samples respectively, which we revisit in this 
paper. All ten of the previously resolved binaries were 
independently recovered with kernel phases, and for all 
we report significantly improved astrometric precision. 
In addition to confirming the technique and software on 
an unambiguous sample, the dramatic improvements to 
the binary parameters offer the chance to determine or- 
bital elements and therefore dynamical masses. Both 
the detections and the remaining unresolved binaries are 
used in quantitative exploration of the performance lim- 
itations of kernel phase analysis in the recovery of high 
contrast systems. 

2.2. Kernel Phase Analysis 

Kernel phase analysis follows the principles introduced 
in Martinache (2010). The first step is to generate a 
model of the pupil of the imaging system as seen from 
the detector. This task is straightforward based on in- 
formation contained in th e Tiny Tim v. 7.2 PSF sim- 
ulation package for NICl (jKrist et al.l l201 ih . available 
at tinytim. stsci . edu. For kernel phase analysis, the 
model pupil is discretized into a square grid array of 72 
sub-apertures with a unit spacing 1/12 th that of the 
pupil diameter (cf. Figure Regions of the primary 
blocked by spiders or the secondary mirror are not sam- 
pled, and one can also observe that the unit baseline 
imposed by this sampling of the pupil imposes an outer 
working angle of 6X/D. 
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This geometry fihs the (li, 'u)-plane with a regular grid 
of 176 distinct sample points at a cadence of 12 points 
across the diameter. The transfer matrix that relates 
instrumental phase errors to spurious ('U,v)-phase in- 
formation is therefore a 72 x 176 rectangular matrix, 
whose SVD reveals 36 non-zero singular values (that is 
exactly one half of the entire number of sample points 
in the pupil), leaving 176 — 36 = 140 independent kernel 
phase relations. Kernel phase analysis can therefore re- 
cover 140/176 = 80% of the available phase information 
present in the quantized grid. 

The specific discretization chosen was found to have 
relatively little impact on the performance of the algo- 
rithm. If instead we adopt a finer-sampled pupil model 
with 20 points across the diameter, we get 1516 kernel 
phases out of 1632 distinct baselines, leading to a 93% 
phase recovery. We analyzed a portion of our dataset 
with this finer 20 point sampling and found little im- 
provement in the quality of fit or precision in parame- 
ter estimation, though we note somewhat better agree- 
ment between H and J bands with the finer pupil model 
using Levenberg-Marquardt model fitting. Because the 
finer grid analysis was computationally expensive but 
yielded only a small change in fitted binary parameters, 
it was judged that application over a large grid to fit the 
available data was unwarranted. We have therefore cho- 
sen the coarser model for our fitting routines, but note 
that more computer time may produce some improve- 
ments with a finer pupil model. For application to wider 
separation binaries, however, the finer model would be 
strictly required: if we have p points across the pupil, 
Nyquist's sampling theorem imposes an outer working 
angle p\/2D. If this condition is not met, the Fourier 
plane fringes will not be well-sampled and parameter es- 
timates will be subject to aliasing or may not be recov- 
ered at all. 

A shorter wavelength of observation delivers an in- 
crease of angular resolution, but with the same level of 
optical aberration, this also precipitates a greater degree 
of image degradation (lower Strehl ratio). When consid- 
ering residual phase noise, we therefore expect that the 
kernel phase signal-to-noise will be accordingly higher for 
images taken at longer wavelengths. 

We are furthermore limited by the fact that we only 
have single snapshots of each target: without multiple 
frames it is difficult to calibrate systematic errors and ex- 
plore statistical uncertainties on the kernel phase observ- 
ables. We therefore selected a sample of stars for which 
we could see no PSF abnormality or obvious Fourier 
phase structure, and repeatedly applied a Levenberg- 
Marquardt fitting algorithm to the raw kernel phase data 
to attempt to find binary companions. Those targets for 
which no companion model was significantly preferred 
over a single source were deemed to be "unresolved" . We 
then used this unresolved population to establish uncer- 
tainties as ensemble standard deviations for each kernel 
phase, which in turn enables quantification of significance 
in subsequent explorations entailing ^ fitting. The re- 
sults presented here could be considerably improved with 
the design of an observational campaign at the outset 
which delivers better diversity, by exploiting multiple ex- 
posures and dedicated point source calibrators. A more 
comprehensive understanding of systematic errors and 
noise estimates for individual targets would yield more 



sensitive Hmits on detection and better errors on fitted 
parameters. 

2.3. Bayesian Methods 

A binary system at any one epoch can be characterized 
by its angular separation ^, position angle and contrast 
ratio c. The likelihood of a binary model with these pa- 
rameters given the set of kernel phases {K(j)j} is related 
to the statistic by 

L((5,^,c|{i^(/),})aexp(-xV2). (1) 

When normalized, this likelihood is the joint density 
probability function for all three parameters. When cal- 
culating we found it necessary to add an additional 
systematic error term in quadrature to bring the mini- 
mum reduced down to 1 . Confidence intervals for any 
individual parameter can be calculated by integrating 
over the two other parameters. After this marginaliza- 
tion, we estimate the parameter and its uncertainty from 
the mean and the standard deviation of the 1-D marginal 
distribution respectively. 

The approach closely follows established practice with 
closure-phase in NRM-interfero metry for th e character- 
ization of binaries (Ma rtinache et al.ll2009f ). When ap- 
plied to our sample HST imaging data set, the final re- 
sults from our algorithm were: (1), confirmation of bina- 
ries already identified with other methods; (2), the deter- 
mination of statistically sound constraints on the binary 
parameters; and (3), a robust statistical estimate for the 
probability that signals extracted from any given system 
betray the presence of a companion or can be attributed 
to noise. 

The sampling and integration of the likelihood func- 
tion given in Equation [1] is in general difficult, and is 
typically performed by a computationally-expensive grid 
integration or a Monte Carlo Markov Chain (MCMC) 
random sampling method. In this paper we apply 
a recently developed alternative, namely nested sam- 
pling. This method, proposed by Skilling (2004), uses 
an unusual change of variables to calculate the model 
evidence. It has recently seen a s urge of interest; 
e.g. for cosmological model fitting (|Mukherjee et al.l 
120061 : iMukheriee fc ParkinsonI [2QQ8D, and the analysis 
of simulated grav itational wave data (|Avlott et al.ll2009l : 
iFeroz et al.ll2009f ). 

The key idea of nested sampling is to populate the al- 
lowed prior space with a large number (~ 100) of 'ac- 
tive points' which are initially chosen at random and 
subsequently evolved towards ensemble states of succes- 
sively higher likelihood usin g MCMC metho d s. Ou r im- 
plementation was based on ISivia fc Skillind (|2006[ ) and 
ultimately yielded a statistical representation of the like- 
lihood space which could be used for binary hypoth- 
esis testing and estimation of model parameter values 
and their uncertainties. Although a number of alter- 
nate gradient-descent and MCMC methods were bench- 
marked, nested sampling was found to be computation- 
ally the most efficient. A global binarity analysis of the 
entire sample in both the J and H bands could be ac- 
complished quickly, however for objects which are in the 
barely resolved limit, there are well known strong pa- 
rameter degeneracies - particularly between separation 
and brightness of a companion. This ambiguity conffates 
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Figure 1. Diagrams showing (left) the discretized pupil model 
used for kernel phase analysis and (right) the resultant (li, v) sam- 
pling points. 

bright close companion models with somewhat more dis- 
tant fainter companion models, considerably diminishing 
the astrophysical utility of the findings. 

In addition to separately fitting image data in J and H 
bands, nested sampling was fast enough to enable joint 
four-parameter fitting of both images simultaneously 
M{S,6,ch^cj). The ambiguity in separation/contrast 
from separate fitting was found to be greatly ameliorated 
by the covariance of separation with contrast. Joint fit- 
ting enforces identical separation between bands, greatly 
restricting the size of the valley of degeneracy with 
contrast ratio. These findings are promising for the com- 
ing generation of Integral Field Unit cameras working 
with extreme AO systems which naturally deliver spatio- 
spectral data cube observations. 

For the joint fitting, an additional error term was 
added in quadrature to represent unknown noise sources. 
This was found iteratively such that each band separately 
had a minimum reduced of 1 at the best joint fit pa- 
rameters. In cases where the existing error estimates 
resulted in a minimum reduced < 1^ iio adjustment 
was made. 

For this study, we searched a delimited parameter 
space for companions. We initially searched up to a 
contrast ratio of 200, somewhat beyond the limits es- 
tablished in Section 13.41 Candidate binaries were then 
compared to Plots [2] and [3] to establish significance. The 
range of separation explored ran from 30 to 200 mas and 
all position angles were considered. 

3. RESULTS AND DISCUSSION 

Images for all objects in Tabled] in the two filter bands 
(FllOW and F170M) were recovered in digital form from 
the HST MAST Archive, where they are listed under 
Proposals 10143 and 10879. All data were processed us- 
ing our kernel phase techniques, and in the discussion 
which follows, we divide our results into three subsec- 
tions: 13. H binaries already reported; 13. 2[ new binaries 
uncovered by kernel phase: [331 marginal detections mer- 
iting further study; and l3.4[ sample detection thresholds 
and the incidence of unresolved sources. 

3.1. Known binaries 

In all c ases where companions were reported by 
iReid et al.l p006a) or Reid et al. (2008), our new anal- 
ysis independently recovered strong systematic signals 
confirming binarity. We therefore confirm all previously 
reported detections, and we stress that our analysis was 
blind in the sense that no prior knowledge was employed 
in our search. 



Table 1 

Sample of u nresolved L - dwar fs from IReid et al.l (|2006al ) and 
IReid et al.l (l2008l ) (after line break) 



2 MASS name 


Sp. type 


J 


H 


K 


2MASS 


J00361617+1821104 


L3.5 


12.47 


11.59 


11.06 


2MASS 


J00452143+1634446 


LO 


13.06 


12.06 


11.37 


2MASS 


J01075242+0041563 


L8 


15.82 


14.51 


13.71 


2MASS 


J01235905-4240073 


M8 


13.15 


12.47 


12.04 


2MASS 


J01550354+0950003 


L5 


14.82 


13.76 


13.14 


2MASS 


J02132880+4444453 


LI. 5 


13.51 


12.77 


12.24 


2MASS 


J03140344+1603056 


LO 


12.53 


11.82 


11.24 


2MASS 


J03552337+1133437 


L6 


14.05 


12.53 


11.53 


2MASS 


J04390101-2353083 


L6.5 


14.41 


13.37 


12.81 


2MASS 


J04455387-3048204 


L2 


13.41 


12.57 


11.98 


2MASS 


J05002100+0330501 


L4 


13.67 


12.68 


12.06 


2MASS 


J05233822-1403022 


L2.5 


13.12 


12.22 


11.63 


2MASS 


J06244595-4521548 


L5 


14.48 


13.34 


12.60 


2MASS 


J06523073+4710348 


L4.5 


13.55 


12.37 


11.69 


2MASS 


J08251968+2115521 


L7.5 


15.12 


13.79 


13.05 


2MASS 


J08354256-0819237 


L5 


13.15 


11.95 


11.16 


2MASS 


J08472872-1532372 


L2 


13.52 


12.63 


12.05 


2MASS 


J09083803+5032088 


L7 


14.56 


13.47 


12.92 


2MASS 


J09111297+7401081 


LO 


12.92 


12.20 


11.75 


2MASS 


J09211410-2104446 


L2 


12.78 


12.15 


11.69 


2MASS 


J10452400-0149576 


LI 


13.13 


12.37 


11.81 


2MASS 


J10484281+0111580 


LI 


12.92 


12.14 


11.62 


2MASS 


J10511900+5613086 


L2 


13.24 


12.42 


11.90 


2MASS 


J11040127+1959217 


L4 


14.46 


13.48 


12.98 


2MASS 


J11083081+6830169 


L0.5 


13.14 


12.23 


11.60 


2MASS 


J12130336-0432437 


L5 


14.67 


13.68 


13.00 


2MASS 


J12212770+0257198 


LO 


13.17 


12.41 


11.95 


2MASS 


J14283132+5923354 


L5 


14.78 


13.88 


13.27 


2MASS 


J14482563+1031590 


L5 


14.56 


13.43 


12.68 


2MASS 


J15074769-1627386 


L5 


12.82 


11.90 


11.30 


2MASS 


J15394189-05200428 


L3.5 


13.92 


13.06 


12.58 


2MASS 


J15525906+2948485 


LI 


13.48 


12.61 


12.03 


2MASS 


J16580380+7027015 


LI 


13.31 


12.54 


11.92 


2MASS 


J17054834-0516462 


L0.5 


13.31 


12.54 


12.03 


2MASS 


J17312974+2721233 


LO 


12.09 


11.39 


10.91 


2MASS 


J17534518-6559559 


L4 


14.10 


13.11 


12.42 


2MASS 


J18071593+5015316 


LI. 5 


12.96 


12.15 


11.61 


2MASS 


J19360262-5502367 


L4 


14.49 


13.63 


13.05 


2MASS 


J20575409-0252302 


LI. 5 


13.12 


12.27 


11.75 


2MASS 


J21041491-1037369 


L2.5 


13.84 


12.96 


12.36 


2MASS 


J22244381-0158521 


L4.5 


14.05 


12.80 


12.01 


2MASS 


J23254530+4251488 


L7.0 


15.51 


14.46 


13.81 


2MASS 


J23515044-2537367 


L0.5 


12.46 


11.73 


11.29 


2MASS 


J002424. 6-015819 


M9.5 


11.86 


11.12 


10.58 


2MASS 


J010921. 7+294925 


M9.5 


12.91 


12.16 


11.68 


2MASS 


J022842. 4+163933 


LO 


13.17 


12.33 


11.82 


2MASS 


J025114. 8-035245 


L3 


13.08 


12.26 


11.65 


2MASS 


J025503. 5-470050 


L8 


13.23 


12.19 


11.53 


2MASS 


J031854. 0-342129 


L7 


15.53 


14.31 


13.48 


2MASS 


J044337.6+000205 


M9 


12.52 


11.80 


11.17 


2MASS 


J083008. 3+482848 


L8 


15.44 


14.34 


13.68 


2MASS 


J085925. 4-194926 


L7 


15.51 


14.44 


13.73 


2MASS 


J102248. 2+582545 


LI 


13.50 


12.64 


12.16 


2MASS 


J102552. 3+321235 


L7 


15.91: 


15.59: 


15.07 


2MASS 


J104307.5+222523 


L8 


15.95 


14.75 


13.99 


2MASS 


J105847.8-154817 


L3 


14.18 


13.24 


12.51 


2MASS 


J115539.5-372735 


L2 


12.81 


12.04 


11.46 


2MASS 


J120358.1+001550 


L4 


14.01 


13.06 


12.48 


2MASS 


J130042. 5+191235 


LI 


12.71 


12.07 


11.61 


2MASS 


J142131. 5+182741 


LO 


13.23 


12.43 


11.94 


2MASS 


J142528.0-365023 


L3 


13.75 


12.58 


14.49 


2MASS 


J143928.4+192915 


LI 


12.76 


12.04 


11.55 


2MASS 


J150654.4+132106 


L3 


13.41 


12.41 


11.75 


2MASS 


J151500.9+484739 


L6 


14.06 


13.07 


12.57 


2MASS 


J172103.9+334415 


L3 


13.58 


12.92 


12.47 


2MASS 


J200250.7-052152 


L6 


15.32 


14.23 


13.36 


2MASS 


J202820.4+005227 


L3 


14.30 


12.38 


12.79 


2MASS 


J214816.3+400359 


L6.5 


14.15 


12.78 


11.77 


2MASS 


J223732. 5+392239 


M9.5 


13.35 


12.68 


12.15 
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In common with closure phases, non-zero excursions 
in the kernel phases encode information about asymmet- 
ric structures, although the abstracted nature of the ob- 
servable makes for significant challenges in intuitive data 
presentation. One approach to present the way a binary 
signal is encoded upon the kernel phases, and the fin- 
gerprint match of this complex function to the actual 
recorded data, is to simply plot the best-fit model quan- 
tities against the observed data in a correlation diagram. 

Figures [6] and [5] present such a diagram of the bi- 
nary model fit against extracted kernel phase data 
for two illustrative cases: previously-known binaries 
2M 0700+3157 and 2M 0147-4954. The one-to-one cor- 
respondence line is overplotted, delineating the locus of 
perfect fit. 

Kernel phase analysis has yielded greatly improved as- 
trometric precision on most previously known binary sys- 
tems. Whereas previous studies, relying on visual analy- 
sis and PSF subtraction, quoted separations to the near- 
est 10 mas and position angle to ^ 1°, kernel phase 
delivers about one order of magnitude better precision. 
All fitted binary parameters agree, typically to within 
Icr, when separate fits to J and H band kernel phases 
are computed. Table [2] gives final best fit parameters 
from the simultaneous J / H four-parameter fit. In many 
cases, the best fits differ from previously published esti- 
mates, sometimes very substantially (althoug h formal er- 
rors were not quoted in Reid et al. (2006a) or lReid et al.l 
^2008)). 

Counterintuitively, for some of the most readily appar- 
ent binaries (e.g. 2M 0004-4044, 0025+4759, 0915-0422 
and 2152+0937), kernel phase methods proved problem- 
atic and could fail to converge to a good fit. This was 
particularly true for well-separated, low contrast systems 
which are most easily discerned by simple inspection. 
Wide binaries induce phase curvature not well sampled 
by our pupil model, and for these cases, a direct fit to 
the squared visibilities was performed. Visibility data 
were calibrated by dividing by the ensemble means over 
the sample, with dimensionless errors of ~ 0.05 added in 
quadrature. As with the kernel phases, statistical anal- 
ysis was based on nested sampling and the results also 
yielded overwhelming improvements in astrometric pre- 
cision. 

The PSF of the individual target 2M 0004-4044 was 
truncated at the edge of the image, presumably due to 
spacecraft mis-pointing, and therefore did not permit a 
useful kernel phase fit in H band. Nevertheless, a J band 
kernel fit w as found to agree well with the parameters 
published in lReid et al.l ()2006a[ ). and likewise a joint vis- 
ibility fit agreed well in both bands. The results of the 
visibility fit are quoted in Table [2l It is unclear, how- 
ever, what effect the data edge truncation may have on 
our interferometric observables and therefore the values 
given in Table [2] are likely to be subject to an additional 
unknown error for this system. 

Example correlation diagrams and a corresponding 
NICMOS image are shown in Figure |4] and Figure [5l 
Note that while the PSFs of the primary and companion 
overlap and are difficult to visually distinguish, they per- 
mit a clear kernel phase fit with very precise parameter 
estimates. 



3.2. Discovery of New Binary Candidates 

Table [3] reports five firm bina r y cand ida tes not de- 
tected in the original Reid et aTl (|2006af ) or iReid et all 
(2008) studies, but recovered at very high 99.9% confi- 
dence from our kernel phase analysis (with the exception 
of 2M 0045+1634 for which the detection confidence was 
only 99%). Correlation plots for all of these are displayed 
as Figures [8] to [TTJ 

Reid et al. (2006a) noted three targets which exhib- 
ited broad PSFs, however they went on to report these 
stars had "... no evidence for the presence of a sec- 
ondary component, and the broader profiles are probably 
an instrumental effect." These were 2M 1507-1627 and 
2M 1936-5502, with a PSF FWHM of 2.47 pixels (106 
mas) and 2M 0036+1821 with a FWHM of 2.56 pixels 
(110 mas), as opposed to the FWHM of 2.3-2.4 pixels 
found through the rest of the unresolved sample. Our 
kernel phase analysis identifies two of these as binary 
candidates: 2M 1936-5502, and 2M 0036+1821. We note 
that for both of these, the kernel phase signal-to-noise ra- 
tio is only of order 2, and they exhibit a correspondingly 
noisy correlation plot. 

2M 1936-5502 supported an alias fit at around 225° 
position angle, and the contrast in H band was very 
poorly constrained. This may well be considered the 
most marginal fit reported in this section. Nevertheless, 
both bands support overlapping position angle modes at 
330°, and this object was considered for further study. 

For 2M 0036+1821, there are two distinct minima in 
H band, one of which overlaps precisely with the single 
distinct minimum from the J band data. In assigning 
parameter estimates for this object, we have assumed 
some uncalibrated source of noise affected the H band 
and have therefore restricted the parameter space de- 
liberately to contain only that region around the high- 
significance J band detection. 

The third target with a reported broad PSF, 2M 1507- 
1627, shows marginally-significant companion fits with 
parameters which are inconsistent between J and H 
bands. Until additional data can be recovered with 
higher signal-to- noise, we classify this object to be un- 
resolved with no companion. 

A companion to 2M 0036+1821 was also reported by 
Bernat et ajj (|2010l ) with non-redundant masking inter- 
ferometry in K band with the laser guide star adaptive 
optics system on the Hale Telescope at Palomar Observa- 
tory. From data taken in September 2008 (some 3.5 years 
after the HST observations), these authors report a con- 
trast ratio of 13.1 ± 3.1 and separation of 90 ± 11 mas 
at a position angle of 114° ± 5. Although changes in the 
binary separation parameters are to be expected with 
progress in an orbit, the markedly different contrast ra- 
tio when compared to our fit in Table [3] is, at first glance, 
hard to reconcile. However we note that both studies em- 
ploy similar fundamental methodologies and are affected 
by the strong separation-contrast ratio degeneracy pr evi- 
ously discussed. In particular, the work of Berna t et al.l 
(|2010[ ) observed only in the Kg band and therefore en- 
joyed none of the advantages offered by dual wavelengths 
in lifting the ambiguity discussed in Section |231 Indeed, 
in their discussion an alternate detection at 25:1 contrast 
and 243 mas with equal probability is debated and ruled 
out based on the HST archival data. To test the hy- 
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Table 2 

Model Parameters for Known L-Dwarf Binaries 



2 MASS 


Sep. 


Pos. angle 


Contrast 


Contrast 


Number 


(mas) 


(degrees) 


Ratio (J) 


Ratio (H) 


0004-4044^ 


84.6±0.2 


224.6±0.1 


1.04±0.01 


1.02+0.02 


0025+4759^ 


329.0±0.3 


233.04i0.06 


1.32±0.04 


1.03+0.04 


0147-4954 


139.8±0.1 


72.66±0.05 


2.37±0.09 


2.06+0.06 


0429-3123 


525.2±1.2 


285.3±0.2 


3.51±0.1 


2.82+0.06 


0700+3157 


179.6±0.5 


105.8±0.1 


4.52±0.07 


3.81+0.03 


0915-0422^ 


738.6±0.15 


26.89±0.01 


1.114+0.002 


1.264+0.002 


1707-0558 


1009.5±1.0 


34.9±0.05 


10.6+0.15 


7.5+0.2 


2152+0937^ 


253.7±0.09 


94.5±0.02 


1.09+0.03 


1.15+0.03 


2252-1730^ 


125.9±0.4 


353.5±0.1 


2.46+0.01 


3.395+0.03 


2255-5713 


178.6±0.4 


172.7±0.1 


5.05+0.08 


4.33+0.02 



^ Low contrast: fit with visibilities. See Section 13.11 
b This object is the subject of .Reid et al.1 (|2006bl 



pothesis that both the iBernat et al.l ()2QlQf ) companion 
and the one reported here are consistent with the same 
degree of phase asymmetry, we re-fit our kernel phase 
data with an enforced higher contrast ratio above 10:1. 
This immediately resulted in much larger bes t-fit separa- 
tions, approaching those from B ernat et al.l (|2Q1 Q), and 
we therefore conclude that the two studies have probably 
identified the same companion. 

All of t hese systems were orig:inally assigned spectral 
classes in iReid et al.l (|2006a[ ). which must be altered to 
refiect the discovery of new companions. A similar ap- 
proach is taken here, usi ng tabulated J — H colours 
(|Pecaut fc Mamajekl lin preparatio n) to determine new 
spectral classes, taking into account the contrast ratio 
found at each wavelength. The class assigned to the 
secondary from its colour is compared to the expected 
contrast with the primary as a function of spectral class, 
and found in each case to be relatively consistent. These 
classes, being based on contrasts which are themselves 
subject to error, are accurate only to within one division. 
Of special interest are 2M 2351-2537 and 2M 1936-5502. 
In the former case. lReid et al.l(|2QQ8 ) reassigned the spec- 
tral type from L0.5 determined in Reid et al. (2006a) to 
M8, and excluded the system from their 20-parsec cata- 
logue of L dwarfs o n these grounds. On the other hand, 
I Andrei et al.l (|2011f ) classify it as LO.O based on precision 
photometry. Revealing it as a binary system, its colours 
and luminosity imply that both primary and secondary 
must be early-L dwarfs, with possible classes LO and LI. 

Given the uncertainty in contrast in the 2M 1936-5502 
system, no accurate spectral type can be determined for 
the secondary component. If the high-contrast fit is con- 
firmed, then the secondary must have a spectral type 
later than T9. A brown dwarf of Y class is possible in 
principle, but coeval companionship to an early-L pri- 
mary seems unlikely if the primary is a brown dwarf. If 
this is the case, then the companion may be of planetary 
mass. On the other hand, if the primary is a star, then 
this could be an evolved companio n of br own dwarf mass 
as discussed i n iKir kpatrick etldl (|2012f ). We note that 
iFahertv et al.l (|2012l) present a new parallax distance 
for this object of 15.08 ±1.2 pc, which con ipares well 
with th e spectroscopic distance estimate in iReid et al] 
(|2006af ) and comparatiely poorly wi th new estimates we 
calculate with t he revised tables of iPecaut fc Mamaiekl 
(|in preparation! ), from which we obtain a distance of 



18.6 ± 2.8 pc, though we note that the measurement is 
still consistent with the quoted uncertainties. 

For 2M 0036+1821, the distance is kno wn from 
trigonometric parallax to be 8.77 ± 0.06 pc (Da hn et al.l 
2002). The angular separation of 44.5 mas therefore 
equates to 0.4 AU projected physical separation. This 
is therefore one of the closest projected separations of 
any resolved brown dwar f binary. The 90 mas separation 
from IBernat et aP (|2010f ) gives typical binary orbital pe- 
riods of ~2-3 years, while on the other hand our 44.5 mas 
separation gives orbits of < 1 year. Ignoring the changes 
in separation (due to degeneracy errors), the position an- 
gle change of 84° observed between the HST and Palomar 
datasets could be consistent with more than one orbit in 
the former case or several orbits in the latter: either is 
consistent with the data and a more rapid observing ca- 
dence is required to unambiguously follow the orbit in 
this particular system. 

Excepting 2M 0036+1821 and 2M 1936-5502, the dis- 
tances of other objects are known from spectroscopic par- 
allax and are therefore subject to change depending on 
whether the system is found to be a binary. In general, 
discovery of a companion implies a higher total lumi- 
nosity for the system, and therefore a greater distance. 
Table [3] lists recalculated spectroscopic distances, with 
standard 15% uncertainties. In particular, the reassign- 
ment of 2M 2351-2537 as an LO/Ll binary moves it back 
into the 20-parsec L dwarf sample, while 2M 0045+1634 
and 2M 2028+0052 fah beyond 20 parsecs. The spectro- 
scopic distance of 9.2 ± 1.4 pc for 2M 0036+1821 agrees 
weh with the spectroscopic parallax of 8.77 pc. 

Our kernel phase re-analysis o f the HST archive d oes 
not detect a binary reported by IBernat et al.l (|2010f ) in 
2M 0355+1133 at 90 mas separation. One explanation 
may be that at the HST epoch the projected separa- 
tion was smaller and therefore unresolved. This is nev- 
ertheless included in the discussion of binarity fraction 
in Sect ion lOl The obj ect 2M 0045+1634 was observed 
by Stumpf et ajj (|2010f ) with differential spectral imag- 
ing using the HST NICMOS camera; a PSF variation 
was noted but no other signal found in the vicinity of 
the brown dwarf. 

Thus four of the five detections reported in Table [3] 
can be considered new, and we confirm the previously- 
reported detection in 2M 0036-hl821 albeit with dis- 
crepant best-fit contrast. All detections are significantly 
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below the formal diffraction limit and at separations 
of the order of the detector plate scale, which is 43.1 
mas/pixel. It is therefore unsurprising that the candidate 
binaries are a subject to parameter correlation and sub- 
stantial uncertainties. These may therefore require con- 
firmation and accurate constraint of parameters found 
with follow-up observations employing a larger telescope. 

3.3. Possible New Companions at Higher Contrasts 

Estimates in 'Martinach^ (j2Qllf ) suggested the possi- 
bility of detections at contrast ratios greater than 50:1 
using kernel phases. The recent di scovery of planetary - 
mass companions to brown dwarfs (iTodorov et al.|[2QTQ[ ) 
indicates that these objects exist and are at least com- 
mon enough that in a sufficiently large dataset candidates 
may be found. On the other hand, given the parameter 
degeneracy, it is possible that some high contrast detec- 
tions revealed by these techniques are in fact very close 
low contrast binaries. Owing to this large uncertainty, no 
attempt has been made to assign separate magnitudes or 
spectral classes to the faint companions identified in this 
section. 

A study of all the unresolved objects in our sam- 
ple extending to a contrast of 100 revealed four ob- 
jects with fits agreeing in position angle between bands: 
2M 0314+1603, 2M 1539-0520, 2M 0830+4828 and 2M 
0109+2949. 

These objects were then subjected to a nested sampling 
joint fit to examine correlation at best fit. In the extreme 
contrast-separation regime probed here, the parameters 
show very pronounced degeneracy and it is impossible to 
reliably distinguish between high contrast objects further 
out and lower contrast objects closer in, even using both 
J and H bands. No error was added in quadrature, as 
in each case the kernel phases were already over-fit with 
reduced < 1- 

2M 0314+1603 is the best of the candidates in Table [5] 
and exhibits excellent correlation in J band. Neverthe- 
less, while H band data favours a minimum at this 
same position angle, it is unconstrained in contrast ratio 
and reliably runs off to high contrast in each attempted 
fit. Likewise 2M 1539-0520 supported good correlation 
diagrams but a surprisingly large discrepancy in contrast 
ratio between bands. This is likely to be the result of a 
strongly degenerate fit which constrains neither contrast 
well; the contrast ratios should therefore not be consid- 
ered well-determined. 

With a J contrast > 160, as a coeval companion to 
a brown dwarf this can only be a planetary mass ob- 
ject. On the other hand, it is possible that the L pri- 
mary may be a star at the low-mass limit of the main 
sequence and its companion is itself an evolved brown 
dwarf. Su ch a system is s till of considerable interest: 
as noted in iDieterich et al.l ()2012f ). low-mass binary sys- 
tems tend to be of equal mass, increasingly so at lower 
primary masses. Accordingly, few brown dwarf secon- 
daries to main sequence primaries are known, especially 
at close separations. This is the low- mass end of the 
'brown dwarf desert' first identified by iMarcv fc Butlerl 
(|2000l ) , which has since become the subject of intense 
study dG rether fc Lineweaver' '2006"; 'Alle n et al.l 12007 : 
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therefore be a potentially significant finding. 

Correlation plots for the marginal detections are shown 
in Figures [12] to [151 

3.4. Survey Confidence and Detection Limits 



HST Co-joint Fit Significance Diagram 




200 300 
Angular separation (mas) 



Figure 2. HST snapshot contrast-detection limits for co-joint fit- 
ting as a function of separation, averaged over position angle. The 
90%, 99% and 99.9% contours are overplotted with labels. Paler 
regions indicate higher significance. 



HST Co-joint Fit Significance Diagram 




Dieterich et al.ii2Q12i : .Evans et al. 2012). A confirmation 
of the status of this system, regardless of result, would 



Angular separation (mas) 

Figure 3. HST snapshot contrast-detection limits for co-joint fit- 
ting as a function of separation, averaged over position angle. The 
90%, 99% and 99.9% contours are overplotted labels. Paler regions 
indicate higher significance. This figure represents a simulation of 
the region near the origin of Figure[2with a finer sample grid. Note 
the turnaround at low contrast: kernel phase performs poorly at 
detecting very-low-contrast companions. 

In order to quantify our survey detection threshold 
as a function of model parameters, we performed a 
Monte-Carlo study simulating detection of a population 
of model binaries. We ran 100 simulations, adding a bi- 
nary at each point on a grid in separation, position angle 
and contrast ratio, in both J and H bands, as a co-joint 
fit. We then added Gaussian noise randomly to each 
of these, with the distribution given by our measured 
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Table 3 

Model Parameters for New L-Dwarf Binary Candidates 



2 MASS 


Spectral 


Distance 


Sep. 


Pos. angle 


Contrast 


Contrast 


J 


H 


Number 


Tvue 


(vc) 


(mas) 


(degrees) 


Ratio (J) 


Ratio (if) 


masf 


maer 


0036+1821^^ 


L3.5 


8.77 ± 0.06 ^ 


44.5±1.2 


198.4±1.3 


1.85±0.3 


2.9±1.1 


12.47 


11.59 


A 


L4 












12.93 


11.91 


B 


L5-6 












13.6 


13.06 


0045+1634^ 


LO 


26.8 ± 4.0 


50.3±0.7 


300.6±2.5 


1.11±0.02 


1.12±0.03 


13.06 


12.06 


A 


LO 












13.75 


12.75 


B 


LO 












13.87 


12.87 


1936-5502^ 


L4 


15.08 ±1.2^ 


67.1±6.4 


330.9±1.0 


17.7±3.9 


35.5±8.1 


14.49 


13.63 


A 


L4 












14.54 


13.66 


B 


(T-Y) 












17.67 


17.53 


2028+0052 


L3 


26.1 ± 3.9 


45.8±1.2 


107.7±1.1 


1.52±0.1 


3.1± 0.5 


14.3 


12.38 


A 


L3 












14.85 


12.7 


B 


L4 










15.3 


13.9 




2351-2537 


L0.5 


17.8 ± 2.7 


63.3±0.3 


348.8±0.3 


2.4±0.1 


2.24±0.2 


12.46 


11.73 


A 


LO 












12.84 


12.13 


B 


LI 












13.79 


13.006 



^ iRe id et aTl (|2006a') listed these objects as having a significant PSF abnormality. A third object 
showi ng such a n abnormality, 2M 1507-1627, does not show a kernel phase binary fit. 
^ .Bernat et al.l (^2010. ) detected this binary with Palomar adapt ive optics aperture m asking. 
^ This object was observed with spectral differential imaging bv lStumpf et all ([20TqI 1. A PSF broad- 
ening was noticed but no other unambiguous signal was detected. 

^ This distance is known from trigonometric parallax to be 8.77 =b 0.06 pc (iDahn et al.l 1200211. 

^ Using the newer spectral class - absolute magnitude tables from iPecaut Mamaiekl 
(|in preparation! ), we revise the Reid spectroscopic distance from 15.4 pc to 18.6 ± 2.8 pc. 



Table 4 

Mass Estimates for New L-Dwarf Binary 
Candidates 



2 MASS 


^0.5Gyr 


MiGyr 


M^Gyr 


Number 


M0 


Mq 


M0 


0036+1821 A 


0.049 


0.064 


0.073 


B 


0.045 


0.060 


0.071 


0045+1634 A 


0.07 


0.078 


0.081 


B 


0.07 


0.078 


0.081 


1936-5502 A 
B 

2028+0052 A 


0.049 


0.064 


0.073 


0.052 


0.066 


0.074 


B 


0.049 


0.064 


0.073 


2351-2537 A 


0.07 


0.078 


0.081 


B 


0.059 


0.072 


0.075 



Table 5 

Model Parameters for Marginal High Contrast Companions 



2 MASS 


Separation 


Position angle 


Contrast 


Contrast 


Number 


(mas) 


(degrees) 


Ratio (J) 


Ratio (if) 


0109+2949 


49 ± 13 


268 + 3 


40 + 25 


43 + 29 


0314+1603 


124 ± 16 


227 + 2 


70 + 10 


> 160 


0830+4828 


48 + 9 


120 + 3 


29 + 16 


14 + 11 


1539-0520 


35 + 5 


332 + 5 


(4 + 3) 


(28 + 16) 



error distribution as discussed in Section 12.21 The de- 
tection rates as a function of separation and contrast are 
averaged over position angle (distributions were found 
to be azimuthally symmetric). Quantitative detection 
thresholds were formulated from comparison of the bi- 
nary model to the null hypothesis. T his method is 
similar to that used in iMartinachd ()2Q1Q| ) . We also ex- 
amined the low-contrast, low separation domain in which 
we have found many of our binaries with a finer sample 



grid. The contrast detection limit curves turn around at 
very low contrasts, as kernel phase performs poorly in 
discerning very- low-contrast companions, which give rise 
to highly-symmetrical images with a weak Fourier phase 
signal. This is true of any phase-based method, as the 
Fourier phases encode information about spatial asym- 
metries in a source image, and this problem therefore af- 
fects closure phases in NRM interferometry as well. This 
therefore places a lower limit on the contrasts detectable 
at very small separations. The contrast thresholds ob- 
tained from these studies are shown in Figures [2] and 
[3] respectively. Our estimates here compare favourably 
with the kernel phase contrast detection limits predicted 
for the Keck Telescopes in Martinachel ((2011), which are 
at best closer to 50:1. This very high sensitivity with the 
HST is unsurprising, given the high quality of the wave- 
front, and this technique therefore holds great promise 
in its application to similar snapshot samples. 

Over any companion detection hangs the question as 
to whether the pair of stars imaged are physically as- 
sociated, or whether they merely happen to lie along 
the same line of sight but are otherwise unrelated. It 
is therefore important to establish the expected count of 
background stars in the direction of each star, which will 
in general differ according to the star's position, owing 
to shape of the Galaxy and the distribu tion of dust. 

The Galaxia software package (Sharm a"et al.l [2QTll) is 
a synthetic survey tool, which calculates the expected 
density of Galactic stars at a given magnitude in a given 
band along a given line of sight. This was applied in 0.1- 
square-degree regions around the coordinates of each bi- 
nary candidate, searching for background stars between 
12th and 18th apparent magnitude in J and H bands. 
For a canonical 12th magnitude primary, these span the 
contrast range from 1 to 100, and therefore cover both 
the brown dwarf and planetary-mass companion regimes. 
This was then scaled down to the expected counts in a 
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200 mas circle around the target. 

In summary, no field in either band shows significantly 
more than 2% probability of finding a background star 
within 200 mas, and in the overwhelming majority of 
cases, including almost all new binary candidates and 
marginal detections, this figure was an order of magni- 
tude lower. It is therefore exceedingly likely that back- 
ground stars do not contribute in any way to the popu- 
lation of companion candidates in this survey. 

4. DISCUSSION 

4.1. Kernel Phase Performance 

From our retrieval of existing and new binaries, we 
have demonstrated that kernel phase interferometry per- 
forms well in medium and wide band filters. For shorter 
wavelengths the PSF quality degrades increasing the er- 
rors, however in all cases space telescope data are firmly 
in the regime where wavefront quality is excellent and 
easily sufficient for the purposes of the algorithm. This 
method has demonstrated the delivery of very precise 
astrometry for medium-separation systems. If employed 
over multiple epochs, this will permit correspondingly ac- 
curate dynamical mass measurements for most detected 
binaries. 

For the very closest companions, we note that the 
signal-to- noise ratio is lower and accordingly kernel phase 
yields weaker constraints on the binary parameters, 
which nevertheless should still permit the determination 
of dynamical masses. For systems whose close separation 
puts them beyond the diffraction limit, errors are domi- 
nated by covariance between the separation and contrast 
creating model ambiguity between close, bright compan- 
ions and distant, fainter companions. We have shown 
that when multi-wavelength observations are employed, 
this degeneracy can be partially lifted. 

4.2. Opportunities for Dynamical Mass and Radius 
Measurement 

Adop ting distance and mass estimates from lReid et al.l 
([2006a) and Reid et al." (2008) , and our own parameter 
estimates, typical binary objects reported here such as 
2M 2351-2537 and 2M 2028+0052 would have a binary 
separation of ~ 1 AU and therefore orbital periods of 3 
to 4 years. A follow-up campaign with LGS AO could be 
used to track the orbits of these binaries and dynamically 
determine their mass in the near term, as has already 
been achieved for GJ 802B in [Ireland et all (I2008D. 

With the alrea dy-observed e pochs from iReid et al.l 
(|2006al . 120081 ) and lBernat et al.l (|2010D . two more epochs 
should be sufficient to permit a ffrst fft to the binary 
parameters of each system. Targets observed in 2006 or 
2008 will have completed a substantial fraction of their 
orbit since they were first observed; on the other hand, 
over a two month period they will have rotated by 6 de- 
grees in mean anomaly, which is substantially greater 
than the ~ degree uncertainties in position angle ob- 
tained with kernel phase or aperture masking. These 
are pessimistic figures, in that the systems are early L 
dwarfs and include several systems substantially closer 
than 2 AU, which leads to correspondingly shorter or- 
bital periods for targets of interest. 

In Table HI we present approximat e mass es com- 
puted using the methods of , Reid et all (|2006a[ ). These 



were calculated by taking the absol ute J magnitude of 
the new spectr al class as listed in iPecaut fc MamajeS 
(|in preparationf) , applying J ma g nitud e bolometric cor- 
rections from iGolimowski et all (j2004D and comparmg 
bolometric magnitude s with the 0.5 , 1 and 5 Gyr 
isochrones from C habrier et al.l (|2000f ). The calculated 
values are subject to direct uncertainties of order ~ 
O.O3M0: despite the significant uncertainty in photom- 
etry, bolometric magnitude is a very strong function of 
mass and age and purely photometric uncertainties are 
small. On the other hand, the models used come with 
significant caveats, especially because brown dwarfs cool 
as they age: the mass estimates, given photometry, there- 
fore depend on the assumed age of the system. Obtain- 
ing dynamical masses is key to calibrating these models 
brown dwarfs in general, as discussed in Section [TJ These 
close new binaries therefore present a significant addition 
to the population of targets which can be followed up on 
short time scales. 

In addition to this, many brown dwarfs are known 
to exhibit periodic radio or Ha emission, modulated 
by the body's ro t ation period, catalogued in detail 
in i Antonova et al.l ( 120131 ). The radio emission is be- 
lieved to be from electron cyclotron maser instability in 
the brown dwarf mag netosphere (^Hallinan et al.l 120081 : 
iKuznetsov et al.l l2012f ). Knowing this rotation period 
and the projected rotational velocity vs'mi from spec- 
troscopic observations, it is possible then to determine 
the radius of the object if its spin is assumed to be per- 
pendicular t o the orbital plane of the system. This has 
been used bv lBerger et al.l (|2009f ) to determine the radius 
of a radio- and H^- variable component of the brown 
dwarf binary system 2MASSW J0746425+200032. This 
new sample of binary systems therefore also presents 
the opportunity to systematically study the radii of any 
brown dwarfs found to have radio emission. Notably, 2M 
0036+1821 has a rotation period of 3.08 ± 0.05 hr as de- 
termined in lHallinan et al.l (|2008f ), and therefore presents 
the clearest new opportunity for radius measurement in 
this dataset. 

4.3. Binarity Fraction in L Dwarf Sample 

We now calculate the revised L-d warf binary f ractio n 
in the 20-parsec sample. Following iReid et al.l (|2008[ ), 
correcting for Malmquist bias requires that for statistics 
in a 20-parsec sample, we should consider only systems 
within 19 parsecs. The unbiased estimator is then simply 
the binary fraction observed. For uncertainties, ordinar- 
ily Poisson statistics would be used; in this case, how- 
ever, the sample is too small (N < 100)^ and so binomial 
statistics are require dJ|Burgasser et al.ll2003f ). 

Where IReid et al.l (p08) had 8 binaries out of 64 sys- 
tems observed in the nearest 20 parsecs, giving = 
12.5^3 0%, five of the previously unresolved systems in 
the sample now support binary detections. Of these, 
the spectral class of 2M 2351-2537 is reclassified to 
LO/Ll, while the increased spectroscopic distance to 
2M 0045+1634 leaves it outside of 20 parsecs. Given 
that 2M 2028+0052 was already known to lie beyond 20 
pc, 3 of these are admissible in the 20-parsec sample. 
Furthermore, 2M 1936-5502 supports a high-contrast 
fit, indicating a potential sub-L-dwarf classification, so 
among the new discoveries, there are 2 firm detections: 
2M 0036+1821 and 2M 2351-2537. In addition to these, 
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we admit a dete c tion o f a companion to 2M 0355+1133 
in iBernat et al.l ()2010l ). for a total of three new bina- 
ries. This yields a new binary fraction of 65 = 17.2^.1'^%. 
If 2M 1936-5502 is included, the binary fraction is then 
lS.75tlr%. 

4.4. Formation of Brown Dwarfs 

A typical brown dwarf has a mass significantly below 
the Jeans mass of a collapsing protostellar gas cloud 
(|Reid fc HawlevI l2005f ), and it is therefore difficult to 
explain the observed abundance of field brown dwarfs 
based on the standard Jeans collapse theory of star for- 
mation. Th is leads to two widely-discussed hypotheses 
(|Basu|[2QT2l ). The first, namely gravoturbulent collapse, 
posits that turbulent gas dynamics allows for the col- 
lapse of smaller clouds than naive analysis suggests, and 
therefore brown dwarfs may form in essentially the same 
manner as other stars. The alternative is embryo ejec- 
tion, whereby gravitational interactions in a protostellar 
system may eject low-mass companion embryos before 
they accr ete enough gas to achieve fusion in their cores 
(|Bate et al. 2002). Bate & Bonneh (2005) propose that 
protostellar cores start out at the opacity limit for tur- 
bulent fragmentation (of order a few Mj) and proceed to 
accrete gas until they are ejected dynamically from the 
cloud. This behavior differs particularly strongly from 
direct turbulent collapse in the substellar regime, so that 
brown dwarfs are thrown out very early in the accretion 
process. Therefore by testing model predictions for these 
ultracool stars we provide one of t he most str i ngent tests 
of the gene ral validity of the Bate fc Bonnelll (|2005f ) and 
iBatd (|2012[ ) models of star formation for all stellar types. 

The chief objection to direct gravitational collapse 
models had previously been that the collapse of such a 
low-mass cloud would require it to start out very dense 
and ve ry cold com pared to well-studied star forming re- 
gions (|Reid fc Ha wlev 2005]^ Recent millimeter inter- 

ferometry by I Andre et al.l ( 2012f ) has caught just such 
a gravoturbulent collapse in the act, observing a grav- 
itationally bound cloud of mass ^ 0.02 — O.O3M0. In 
addition this, iMonin et al.l (|2013f ) report the detection 
of a molecular outflow from the brown dwarf binary FU 
Tau, which is the third young brown dwarf system found 
undergoing a formation process analogous to low-mass 
stars. With at least one example of a brown dwarf form- 
ing directly, models that do not account for this behav- 
ior must therefor e be r efined. Indeed, simulations by 
lJumper fc Fisherl (|2012[ ) reproduce the distribution of 
wide brown dwarf companions to main sequence stars 
as well as the population of close binaries. 

On the other hand, while embryo ejection can easily 
account for the formation of isolated field brown dwarfs, 
binaries are problematic. The binarity fraction of very 
low mass stars can discrim inate between these two mod- 
els (|Burgasser et al.l [20071 ). Embryo ejection predicts a 
low binarity fraction in field brown dwarfs, as the grav- 
itational interaction required would disrupt all but the 
most tightly bound binaries. The predicted binarity frac- 
tion is thus < 5% in t he oldest models for embryo ejec- 
tion ([Bate et al.ll2002f ), rising only to 8 ± 5% for systems 
with a primary mass in the range of 0.08 to 0.1 Mq in 
the most recent models incorporating radiative feedback 
([Bate. 20 12.) . The observed binarity fraction of ultracool 



stars seems closer to 15% (|Bouv et al.l 120031 : iReid et all 
2006a, 2008J ) whic h puts it into confi i ct wit h the models o f 
iBate et al.l (|2002l ): iBil^fc Bonnelll (|2005h : (I20T2I) . 
although the authors noted that this prediction is subject 
to significant uncertainty. 

The high binarity fraction reported in this paper is still 
further evidence in favour of a higher binarity rate than 
predicted by embryo ejection. Moreover, with the addi- 
tional detections this prediction now lies several a away 
from the observed value. The hybrid ejection model of 
iBasu fc Vorobvovl ([2012), by which still-collapsing clouds 
can be ejected from a protostellar disk, makes few clear 
predictions regarding binarity, other than that there 
sh ould be few la rge-separation binaries. Given the results 
of I Andre et al.l (^2012.) in observing a gravitationally- 
collapsing brown dwarf mass object, it is not implausible 
to suggest that direct gravitational collapse may account 
for most or all field brown dwarfs. 

The effect of radiative feedback in amelioratin g the 
difficulties of modelling brown dwarf formation in iBatd 
(j2012D has been suggested by iDieterich et~all (j2012D as 
an explanation for why the initial mass function and com- 
panion mass function changes at the hydrogen burning 
limit; while this idea is highly speculative, it seems likely 
that the identification of systems with a late main se- 
quence pr imar y and a brown dwarf secondary as reported 
in Section [X3l will be of interest in testing this hypothesis. 

5. CONCLUSIONS 

We have shown that extraction and fitting to the self- 
calibrating kernel phase interferometric observables al- 
lows for a significant increase in the robust detection of 
companions and in the accuracy of best-fit parameters re- 
covered when compared to simple inspection of images. 
The Hubble Space Telescope^ or any instrument deliv- 
ering good wavefront quality may benefit dramatically 
from such an approach. 

Using the kernel phase technique on all 79 ultra- 
cool dwarfs in the combined HST samples of iReid et al.l 
(2006a) and Reid et al. (2008), we independently recover 
all 10 prior detections and improve on the precision of fit- 
ted parameters such as position angle and separation by 
a factor of ~10. Furthermore, we report 5 new binary 
detections missed by the original authors, 4 of which are 
presented here for the first time, and additionally four 
marginal detections of close or high contrast compan- 
ions. This population forms an excellent base for dynam- 
ical studies to establish masses; a prospect particularly 
favored by the improved precision in parameter estima- 
tion. As well as these confident new detections, kernel 
phase identifies up to four more marginally-resolved close 
or faint companions, which may be of planetary mass. 

The finding of a larger binary fraction for this sam- 
ple helps shed light on the formation mechanisms of very 
low mass objects. The five additional ultracool binaries, 
if confirmed, lend further support to gravoturbulent col- 
lapse models for the formation of low-mass stars in the 
field. 
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Figure 4. Log scale image of 2M 0147-4954 in F170M filter. The corresponding correlation diagrams are displayed in Figure O Note that 
while the PSFs are hard to distinguish visually, the kernel phase fit is excellent. 
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Kernel Phase Correlation Diagram, target 0147-4954 , J band Kernel Phase Correlation Diagram, target 0147-4954 , H band 
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Figure 5. Correlation diagrams for the previously-detected binary 2M 0147-4954 in (left) J band and (right) H band. 
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Figure 8. Kernel phase correlation diagrams for the newly- detected binary 2M 0045+1634 in (left) J band and (right) H band. 
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Figure 9. Kernel phase correlation diagrams for the newly-detected binary 2M 1936-5502 in (left) J band and (right) H band. 
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Figure 10. Correlation diagrams for the newly-detected binary 2M 2028+0052 in (left) J band and (right) H band. 



Brown Dwarfs from Kernel Phases 




Model Kernel Phases Model Kernel Phases 

Figure 11. Correlation diagrams for the newly-detected binary 2M 2351-2537 in (left) J band and (right) H band. 




Kernel Phase Correlation Diagram, 2M 0314+1603, J band Kernel Phase Correlation Diagram, 2M 0314+1603, H band 




Model Kernel Phases Model Kernel Phases 

Figure 13. Kernel phase correlation diagrams for marginal detection 2M 0314+1603 in (left) J band and (right) H band. 



Kernel Phase Correlation Diagram, 2M 1707-0558, H band 
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